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Abstract 

The possibility that the present acceleration of the nniverse is driven by a kind of viscous fluid 
is exploited. At background level this model is similar to the generalized Chaplygin gas model 
(GCGM). But, at perturbative level, the viscous fluid exhibits interesting properties. In particular 
the oscillations in the power spectrum that plagues the GCGM are not present. Possible fundamental 
descriptions for this viscous dark energy are discussed. 
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1 Introduction 

The observations of the dynamics of galaxies, cluster of galaxies and of the supernova type la indicate 
that about 95% of matter content of the universe is not composed of baryons mm A fraction of about 
1/3 of this dark component appears in the agglomerated structures, and it is called dark matter, while 
the remaining 2/3 appears as a smooth component, driving the acceleration of the Universe, and is called 
dark energy. There is a large number of models trying to take into account the presence of the dark 
component. The most popular one is the so-called ACDM mm, where dark energy is represented by 
the cosmological constant, while dark matter is composed of WIMPS, a cold dark matter, composed of 
weakly interacting massive particles which must be relics of a grand unihed phase of the universe, like 
axions. The ACDM has achieved great success in explaining the observational data (even if there claims 
in the opposite sense 00), but faces at same time many theoretical difficulties, like a huge discrepancy 
(of about 120 orders of magnitude) between the predicted and the observed values of the cosmological 
constant [ZlIHl. 

Other models are competitive with ACDM, like quintessence 0[ini and K-essence HU. Another 
quite recent proposal is the Chaplygin gas mini, which is phenomenologically represented by a fluid 
with negative pressure which varies with the inverse of the density. The Chaplygin gas model has been 
generalized by considering that the pressure, besides to be negative, depends on an arbitrary power of the 
inverse of the density. One of the great advantages of the Generalized Chaplygin gas model (GCGM) is 
the possibility of unifying the description of dark energy and dark matter Olini: a fraction of this exotic 
fluid can agglomerate locally, while the other fraction remains a smooth component. However, criticisms 
have been addressed to this proposal mainly due to its behaviour concerning density perturbations, which 
exhibits large oscillations in the resulting power spectrum which do not appear in the observed power 
spectrum of mass agglomeration m 

In our point of view, the question of the oscillations in the matter power spectrum in the GCGM is 
controversial. The oscillations in the power spectrum of the GCG are not transferred to the baryonic 
power spectrum, and after all, the direct observation concerns baryons. But, there are claims that the 
oscillations in the dark component is reflected in the erg normalization m In this sense, it should be 
interesting to find a way out to this problem, keeping at same time the advantages of the GCGM. 

In this work we will explore the possibility that the present acceleration of the universe is due to a kind 
of viscous fluid. It is well known that bulk viscosity can generate an acceleration expansion |17j . But, 
to our knowledge, such possibility has been investigated only in the context of the primordial universe. 
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concerning also the search of non singular models imiini. We will consider a simple bulk viscosity model, 
in the context of the Eckart formalism m- Naturally, this approach is phenomenological. Moreover, 
the Eckart formalism is not completely consistent, being a truncation of a causal theory 1211123 Hence, 
everything that will be developed here must be later complemented by a fundamental model which can 
justify our phenomenological approach. 

At background level, the description is equivalent to the GCGM: the viscous dark energy interpolates 
a matter dominate phase and a cosmological constant phase. Hence, all observational tests that concerns 
the background behaviour (like the supernova type la data) can be immediately used through the analysis 
already made for the GCGM [221 ■ However, at perturbative level, new features appear: the oscillations 
that plagues the GCG are absent here for a large range of the parameters. We make a simplified 
comparison with the 2dFGRS ignoring for the moment the presence of baryons. The goal is to show that 
it is possible qualitatively reproduce the general features of the mass power spectrum, with the absence 
of expressive oscillations. 

This paper is organized as follows. In the next section we describe the viscous model and deter¬ 
mine under which conditions it can account for the dark component of the universe. In section 3, a 
perturbative analysis is made, and the predicted power spectrum is compared with observations. In 
section 4, we present our conclusions, with some perspectives to a more fundamental motivation of this 
phenomenological model. 


2 Background model 


Let us consider a homogeneous and isotropic universe filled by a fluid with bulk viscosity. For simplicity, 
it will be supposed that the geometry is given by the flat Robertson-Walker metric, 

= dt^ — a^{t)[dx^ + dy"^ + dz^] , ( 1 ) 

where a{t) is the scale factor. The equations of motion are 



p + i% + p*) 
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P* = 


SnGp , 

0 , 

p-mp)- 
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We will consider that p = f3p and that the viscosity coefficient behaves as ^{p) 
equations of motion reduce to 
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From equation we obtain the relation 
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so that the equation can be rewritten as 
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The equation admits the solution 


P = 


Co 


B 


1 + /3 1 + /3 


( 10 ) 


where B is an integration constant and r = 2t[v — 1/2)(1 + (3). 

In the case of the generalized Chaplygin gas, where there is no viscosity, the pressure is given by 


A 


P = 


The relation between the density and the scale factor is given by 

1 

B 


P = 


A- 


a3(i+“) 

The GCGM and the viscosity model coincides, at background level, if /3 = 0 and if 


^ = -la + - 


( 11 ) 


( 12 ) 


(13) 


In this case, initially the Universe behaves as in the matter dominated phase, becoming later dominated by 
a cosmological term. In the case of the viscosity model, the initial phase is characterized by a domination 
of a fluid with equation of state p = (3p, when v < 1/2. For v > 1/2 the behaviour sketched above is 
inverted: initially there is a superluminal expansion followed by a subluminal expansion. The analysis 
based on the supernova type la data for the GGGM can be directly transferred to the dark viscous model, 
since it depends on the background only. In m a extensive analysis of the GGG parameters has been 
made. It has been found that the prediction for the parameter a is a = —0.75/1 q' 24. Hence, positive 
values of (’’normal” viscous behaviour) is preferred. But the dispersion is quite large. 


3 Perturbative study 

In principle, the most interesting situation in the model described above, in view of the present accel¬ 
eration of the universe, concerns the choice (3 = 0 and v = —{a + 1/2). These choices lead exactly to 
the same behaviour of the GGGM for the evolution of the background. However, here we have a more 
normal situation, where the viscosity grows with density when 0 < v < 1/2. On the other hand, most 
of the criticism on the GGGM concerns the fluctuations on the power spectrum which leads apparently 
to a CTg normalization that is not consistent with observation. Hence, in order to verify if the viscosity 
model can lead to improvements with respect to the GGGM, the scalar fluctuations must be studied. 


This study will be done here in the synchronous gauge. 

In order to perform this perturbative study, the field equations are rewritten as 

Rtiu = , (14) 

= 0 , (15) 

= {p + p*)u^u’' -p*g^'' , p* =p-ap)u^-,n ■ (16) 

The equations il4ll5lid are perturbed by introducing 

9)1.1^ = 9)1.1^+ , p = p + 5p , , p*=p* + Sp* , (17) 


where p, and p* are the background solutions described before, while Su^ and Sp* are 

small perturbations around them. The synchronous gauge condition h^o = 0 is imposed. A long but 
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straightforward calculation leads to the following perturbed equations: 
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In these equations, the following definitions were used: 


, hkk . Sp 
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0 = art* 
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( 20 ) 

( 21 ) 


p 

To deduce these equations, the expression for the perturbation of the effective pressure has been used. 
From m, the perturbation in the effective pressure (a crucial aspect for the results to be present later), 
is 


Jp* = /3 - 3- pp - ^(p) 0 - - 


( 22 ) 


These equations can be rewritten in terms of the redshift variable z = —1 + 1/a, where the scale factor 
has been made equal to unity today, oq = 1. Performing also a plane wave expansion in the perturbed 
quantities such that 

5k{t)e^^'^ (fk , k = wavenumber of the perturbation , 


we end up with the following perturbed equations: 
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The primes mean derivatives with respect to the redshift z. The following definitions were also employed: 

(27) 

(28) 

(29) 

(30) 


A = 


1 H” /3 


SttG iy-l/2 

^-CoPo 


4 

































Figure 1: Behaviour for A = 0.1 and i/ = —0.5, 0.3, 0.0, 0.3 and 0.4. The ordinate represents log^Q Fk and 
the abscissa log^Q kh~^. As v grows, the theoretical curve approaches the observed curve. 


The parameter fco is associated with the Hublle length, ko = 2 'kHq/c ~ 27r/(3/i) x 10“^Mpc“^, Hq being 
the Hubble’s constant. The recent results from the WMAP measurements of the anisotropy of the cosmic 
microwave background radiation indicates /i ^ 0.7 |24|. 

We allow the perturbed equations to evolve from z = 500 to z = 0, where the final spectrum is 
computed. The initial conditions are fixed by using the transfer function 


T{k) 


BVk 

TTlkTW 


(31) 


where H is the total density fraction, with respect to the critical density, which in the present case is 
H = 1. The amplitude B can also be fixed by using the normalization of the anisotropy of CMB. Following 
we adopt B — {24 h ^ Mpc)^. At z = 0 we compute the power spectrum 

Pk = \Sk? ■ (32) 


The spectrum is computed for a large range of values of k. The comparison with the observational 
results for the power spectrum of mass agglomeration obtained through the 2dFGRS program is displayed 
in figures 1, 2 and 3 for different values of the parameters v and A. The main feature is the absence of 
oscillations in the power spectrum of the viscous dark fluid. The absence of oscillations occurs for positive 
and negative values of the parameter v. In principle this seems to be surprising, since negative values of 

V should correspond to a positive sound velocity, which should drive oscillations in the power spectrum. 
The reason why these oscillations do not appear is that the spatial gradient term, which drives oscillations 
or instabilities depending on its overall sign, is now composed of three terms, containing not only the 
density contrast, but the velocity and the metric perturbations. The presence of this combination of 
terms avoids the appearance of strong oscillations or instabilities. This combination is due to the form 
of the effective pressure. Just to compare, the GCGM, at perturbative level, contains only the density 
contrast in the spatial gradient term. 

The inspection of figures 1, 2 and 3 reveals that the theoretical curves approach the observational 
data when v ^ 1/2 and A becomes small. A quite reasonable agreement is, for example, obtained when 

V = 0.4 and A = 0.1 (figure 1). Even if the fittings displayed in figures 1, 2 and 3 in general do not reveal 
a remarkable agreement between theory and observation, except in the limits stated above, this is not a 
serious problem in present context due to one fundamental reason: we have not considered baryons. In 
fact, what we should compute is the power spectrum for the baryonic component, perhaps with a bias 
factor which may take into account a contribution of a fraction of the viscous dark fluid. But, what 
we would like to stress is that there is no blow up in the perturbations of this viscous dark fluid, as it 
happens with, for example, the GCGM. Notice that there is a significant suppression of power in the 
spectrum for negative values of v, this suppression is much less important for positive v. Such suppression 
may be interesting if we remember that we are computing the power spectrum of the dark component. 
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Figure 2: Behaviour for A = 0.4 and ly = —0.5, 0.3, 0.0, 0.3 and 0.4. The ordinate represents logio Pk and 
the abscissa logio kh~^. As v grows, the theoretical curve approaches the observed curve. 



Figure 3: Behaviour for A = 0.7 and ly = —0.5, 0.3, 0.0, 0.3 and 0.4. The ordinate represents logio Pk and 
the abscissa logio kh~^. As ly grows, the theoretical curve approaches the observed curve. 
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since the dark component does not agglomerate completely. Hence, the suppression of power in the dark 
component may avoid serious discrepancies with the dynamics of clusters of galaxies. 

4 Conclusions 

In this work, we have developed a phenomenological model for dark energy based on a viscous dark fluid. 
The approach is very simplified since we consider the bulk viscosity in the Eckart formalism, ignoring 
consequently problems of causality. It has been showed that, with the hypothesis that the bulk viscosity 
depends on a power of the density, ^ , interpolation between a matter dominated phase and a 

cosmological constant phase is achieved if < 1/2. Hence, such interpolation can be obtained for a non 
exotic viscous fluid where the viscosity decreases with the decreasing of density. Moreover, the behaviour 
characteristic of the Chaplygin gas model is recovered for negative values of v. 

The evolution of density perturbations for this viscous dark fluid has been computed. There is no 
oscillation in the power spectrum, in opposition to what happens with GCGM. The spectrum is highly 
suppressed for negative values of v but reproduces qualitatively the observed power spectrum for mass 
agglomeration for v positive. The fitting of the observational data becomes quite good when v ^ 1/2 
and H > 0. The model studied here contains just one fluid, the viscous dark energy. Hence, we can 
expect that the adding of baryons will allow to fit reasonably the observational data. The absence of 
oscillations is due to the fact that the spatial gradient of the pressure presents a competition between all 
perturbed quantities. This is dictated by the covariant representation of the bulk viscosity. 

The phenomenological approach developed here must of course be supplemented by a fundamental 
description of this viscous fluid. To do this a specific fluid model must be considered, with some interaction 
between the particles composing this fluid. Topological defects (cosmic strings, domain walls, textures) 
can lead to cosmological fluids with negative pressure in the perfect fluid approximation. We can think 
for example on the evolution of domains wall with friction in an expanding universe m However, the 
effective equation of state for these objects becomes more complex if interactions are taken into account. 
Since interactions are inevitable in a gas of topological defects, it can be expected that deviations from 
the simple perfect fluid approximation can appear. We intend to explore this possibility. 
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